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Abstract. - We study a class of cyclic Brownian heat engines in the framework of finite-time 
thermodynamics. For infinitely long cycle times, the engine works at the Carnot efficiency limit 
producing, however, zero power. For the efficiency at maximum power, we find a universal ex- 
pression, different from the endoreversible Curzon-Ahlborn efficiency. Our results are illustrated 
with a simple one-dimensional engine working in and with a time-dependent harmonic potential. 
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] Introduction. — According to Carnot, the second 
4aw leads to a maximal efficiency 77c = 1 — Tc/Th for 
heat engines working between two heat baths at tem- 
peratures Th > Tc- However, this efRciency can only be 
Achieved in the quasistatic limit, where transitions occur 
infinitesimally slowly and hence the power output van- 
ishes. For real lieat engines, it is more meaningful to 
calculate the efficiency at the maximal power output of 
the machine. The Curzon-Ahlborn efficiency at maximum 
power 77CA = 1 — -^Z Tc/Th was originally derived for a 
Carnot-engine with finite thermal resistance coupling to 
the thermal reservoirs [1]. This behaviour is even recov- 
ered if path optimization techniques are used to maximize 
the work output with respect to the driving scheme, see [2] 
dnd references therein. Recent extensions [3,4] of this ap- 
parent universal law to an infinitesimal series of coupled 
heat engines operating in a linear response regime (i.e. be- 
tween heat baths with small temperature difference) have 
again risen the question whether the Curzon-Ahlborn re- 
sult is quite generally a bound on the efficiency of heat 
engines working at maximum power output. 

In contrast to the macroscopic heat engines considered 
in conventional endoreversible thermodynamics, thermal 
fluctuations play a crucial role in most biologically rele- 
vant systems and dynamics cannot be described on a de- 
terministic (macroscopic) level. In this regime, models for 
thermodynamic machines must incorporate fluctuation ef- 
fects and thus allow also for backward steps even in a di- 
rected motion [5,6]. Brownian motors, in contrast to gen- 
uine heat engines, are mostly driven by time-dependent 



potentials or chemical potential differences. In the last 
two decades, a variety of aspects of such Brownian motors 
has been studied, including dynamics [6,7], stochastic en- 
ergetics [8-11] and efficiencies [12-14], both on a continu- 
ous (Langevin equation) and a discrete (Master equation) 
state space. Since biological motors are typically driven 
far out of equilibrium, linear response thermodynamics is 
not appropriate to describe these systems. The seminal 
work of Sekimoto [8, 9] opened the door for a thermo- 
dynamic description of Langevin systems driven far out 
of equilibrium. Thermodynamic quantities such as work, 
heat, internal energy, and entropy can even be defined on 
a single stochastic trajectory [9, 15,16], yielding the re- 
spective ensemble quantities after averaging. 

Brownian heat engines have been investigated with a 
model system using a spatial temperature profile and a 
ratchet potential [17-19]. This thermal ratchet model has 
also been extended to a Brownian particle in the under- 
damped regime [20] . Efficiencies of such ratchet heat en- 
gines [8, 21-23] and related Brownian heat engines [24] 
have been calculated. Recent studies on the efficiency at 
maximum power, however, have either been restricted to 
the linear response regime [23] or use a questionable def- 
inition of work (which allows for the efRciency to exceed 
the Carnot bound) [22]. For a ratchet heat engine on a 
discrete state space, a violation of the Curzon-Ahlborn 
bound has been found [25]. 

Although widely used as a model system for a stochas- 
tic heat engine, ratchet heat engines are hard to realize 
experimentally because of the necessary temperature gra- 
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dients on a small length scale. Moreover, a macroscopic 
limit of such heat engines does not lead to heat engines of 
the Carnot type. In order to overcome these difficulties, 
we introduce a simple model system of a stochastic heat 
engine whose degrees of freedom are subject to a time- 
dependent potential with time-dependent coupling to and 
decoupling from the two heat baths. This model is some- 
what distinct from the thermal ratchets which are driven 
by spatial temperature differences. Ratchet heat engines 
with a time-dependent but spatially constant temperature 
have been studied in [26-28]. 

We first describe the general setup of a stochastic heat 
engine consisting of two adiabatic and two isothermal 
steps. We then illustrate this setup by a simple one- 
dimensinal example and show that our model exhibits a 
universal dependence of work and heat on the cycling time. 
This leads to a universal form of the efficiency at maxi- 
mum power. We first neglect the effect of kinetic energy in 
our overdampcd description and discuss its role (cf. [13]) 
at the end. 

The model. — We consider a thermodynamic ma- 
chine with internal degrees of freedom x and an exter- 
nal time-dependent potential V{x.,t). On a microscopic 
scale, it is crucial to include thermal fluctuations in the 
description and thus to consider the probability density 
p(x, r) to find the system in state x at time r. In the over- 
damped limit, the time evolution of p(x, r) is governed by 
the Fokker-Planck equation 

drp{^, r) = - V • j = - V • M [- (VV) - keTV] p(x, r) (1) 

where /i is the mobility matrix. Similar to a macroscopic 
Carnot engine, we consider the machine cycle to consist 
of two isothermal and two adiabatic transitions. 

Isothermal steps The potential V{x, t) is varied dur- 
ing a time-interval ti < t < tf at a, given constant tem- 
perature T. The mean work W spent in this process (or 
extracted for W < 0) 



rt f dv 

W — / dr dx p{x,t)-^ 
= Wi„ - TAS + AE 



(2) 



becomes a functional of the time dependent potential 
F(x, r). For notational simplicity, we set fc^ = 1 in the 
following. In Eq. (2), we have introduced the (mean) 
irreversible work 



dr / dx 



j(x,T)^ ^j(x,T) 

p(x,t) 

the (mean) internal energy 

Eir) = J d.Tp(x,T)^(x,T), 
and the system entropy 

S{t) = — dx p(x, r) ln(p(x, r)). 



(3) 



(4) 



(5) 





® 











\ / Pa 







i 






® 








Fig. 1: Sclieme for a cyclic process of a stochastic heat engine 
operating between heat baths at temperature Th and Tc. The 
dotted curve shows a one-dimensional potential V{x,t) and 
the filled curve the boundary values Pa{x), Pb{x) of the cor- 
responding distribution p{x,t). Steps 1 and 3 are isothermal 
whereas steps 2 and 4 are adiabatic. 



see also [15, 16]. 

Adiabatic steps .- Adiabatic steps are somewhat more 
difficult to conceive realistically (just as for macroscopic 
heat engines). Here, they are idealized as sudden jumps 
of the potential while uncoupling from one heat bath and 
coupling to the other heat bath occurs instantaneously 
and thus without heat exchange. Note that, within our 
overdamped description, we first neglect kinetic energies. 
Considering the kinetic energy of the degrees of freedom 
would lead to an additional heat flux in these steps because 
of the relaxation of the momentum due to the tempera- 
ture change. This issue is discussed in detail at the end 
of this Letter. Since there is no time for relaxation, the 
distribution p(x, r) and thus the system entropy S does 
not change during these steps and hence these transitions 
are also isentropic. The mean work applied in such a tran- 
sition occuring at time r is given by 

W ^AE^ J dx p(x, r) [y (x, r + 0) - V{x, t - 0)] . 

(6) 

Cyclic process.- We construct a cyclic process by per- 
forming subsequently the following four steps, sec Fig. 1. 

(1) Isothermal transition at temperature Th during time 
interval < r < <i. The potential V{x,t) is changed 
timc-dependcntly and work is extracted {W < 0). 

(2) Adiabatic transition (instantaneously) from temper- 
ature Th to temperature Tc at time r = ti. The distribu- 
tion Pb(x) = p(x, T = ti) does not change during this step. 

(3) Isothermal transition at temperature Tc during time 
interval ti < t < ti +1^. The potential V{x,t) is changed 
timc-dependcntly (with VF > 0). 
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(4) Adiabatic transition from temperature Tc to tem- 
perature Th at time r = ^1+^3. The distribution 
Pa(x) = p(x, r = ti + ta) docs not change during this 
step. 

For a cyclic process, periodicity in time requires p(x, ti + 
^3) — pI^jO) thus imposing a constraint on the admissible 
time-dependence of V{x,t). Moreover, a cyclic process 
implies = ELi ^S^'^ = 0- Since the sys- 

tem entropy S does not change during the adiabatic steps, 
AS'(2) = AS^^^ = 0, we even have AS'(^) = -Afi-^^) = AS*. 
The total work then is given by 



W 



-I- M^. 

iiT ' irr 



1=1 
(3) 



(7) 



The heat uptake during one complete cycle from the hotter 
heat bath at temperature Th is 



(1) 



(8) 



Example. — In order to illustrate this general concept 
of a stochastic heat engine, we consider the motion of a 
colloidal particle in a trap with time-dependent strength. 
The corresponding potential reads 



V{x,t) = \{t)x^/2. 



(9) 



The probability distribution p{x, r) remains Gaussian for 
all times r if it is so initially. We first derive the equation 
of motion for the variance w{t) = (x^(t)) 



It; = — 2/.tAu' + 2fiT 



(10) 



by multiplying eq. (1) with x'^ and integrating over x, 
where wc denote time-derivatives by a dot. For any initial 
Gaussian distribution p[x,0), the internal energy of the 
system is given as 

E{t)= [ dx p{x,t)V{x) =w{t)X{t)/2. (11) 



Using (10), the mean work (2) can be cast in a local func- 
tional of the new variable w and its first derivative by 
solving (10) for A. For a constant temperature T, this 
yields 



W[X{t)] 



drX 



w 



TAS + AE 



(12) 



_ / dT T\\uwf/ + - \\w]'/ . 



We now have to choose a protocol A(t) for the isother- 
mal transitions. The variance w of the particle position de- 
pends on this protocol via eq. (10). As a consequence, the 
mean work (12) will also depend on this protocol. We ask 
for the optimal protocols A* (r) and A3 (t) yielding a min- 
imal total work (equivalent to a maximal work output). 
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Fig. 2: Stochastic heat engine in a harmonic potential: Time 
dependence of the variance of the particle position w, the tem- 
perature T, the trap strength A and the effective temperature 
Tofi (Eq. (33)) for one cycle with ti = tg — 3, Wa = 0.5, wt = 1, 
Th = 2, Tc = 1, and /i = 1, starting with step 1. 



Rather than imposing boundary conditions on the control 
parameter A, wc impose boundary values of the particle 
position variances Wa = w(0) = w{ti + t3) and Wb = w{ti) 
for a given finite cycle time ttot = ti + t^. While it might 
seem more natural to impose a range of (experimentally) 
possible A values, such an approach can lead to unphysical 
optimal processes which are driven infinitely fast with in- 
finitesimally small changes in the variance w{t) (data not 
shown) . With slight modification of the results in [29] due 
to the different boundary conditions, we get the optimal 
protocols for the variances 

=Wa{l + Wwb/wa - l)rAi)' (13) 

:(t) = Wb (1 + Wwalwb - l)(r - ti)A3)' • (14) 



The optimal protocols A*(r) for the strength of the trap 
can then be calculated from Eq. (10). The resulting quan- 
tities are shown in Fig. 2. 

For this cyclic process, we get the total work output 

-w = i(Th-re)in^--f^ + ^Vv^-v^)^ 

2 Wa H \ti tsj 



in - T,)AS - ( ^ + ^ ) ^hr, 



(15) 



with the irreversible "action" Ai^r- The work only depends 
on the "reduced time" U = l/il/h + l/ts) = tits/iti+ts). 

The heat uptake during one complete cycle from the 
hotter heat bath at temperature Th is given as 



Q 



(1) 



ThAS* - W: 



(1) 



(16) 
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yielding the efRciency 
-W 



(17) 



-Q(l) ThAS - Arr Al 

For ti,i3 oo, we recover the Carnot effieiency rjc = 
1 — Tc/Th. This shows that the sudden adiabatic steps do 
not spoil the reversibility. 

We immediately recognize from Eq. (15) that there ex- 
ists a "stall time" 

(18) 



(Th - T,)AS M(7^h - Tc) In (wb/wa) 

where the heat engine does not perform work anymore. 
Rather, at smaller times i,- < tf', it consumes work. 

We next explore how fast the heat engine must be 
driven to yield a maximal power output. The power 
P = —W I (t\ + ta) has a maximum for 

= (T,-1;W ^^^^ 

which yields a reduced optimal time 

K = f\ = 2tf . (20) 

A similar relation between the stall force and the force at 
maximum power has been obtained in [3] within a linear 
response treatment. The total work output at maximum 
power 



-W* 



1 



(Th - Te)A5 



(21) 



(22) 



(23) 



and the corresponding heat uptake 

_Q*(i) ^ lA5(3rh + Te) 

then yield the efRciency at maximum power 
, ^ 2(rh - Tc) ^ r;c 

m + T, 2-77C/2' 

This result is independent of the imposed boundary val- 
ues Wa and Wb, suggesting a certain degree of universal- 
ity confirmed below. For large temperature differences, 
Th - Tc > Tc, the efficiency has the limit r]* 2/3. 
An expansion for small relative temperature differences, 
r = (Th - T,)/T, « 1, 

if = = -- — + — + (r^) (24) 

' 3r + 4 2 8 32 ^ ^ 

shows that the deviation from the Curzon-Ahlborn effi- 
ciency 

10T3 
2 8~ ^ 32 
occurs at the third order in the relative temperature differ- 
ence T. The validity of the Curzon-Ahlborn efficiency in 
a linear response regime has been shown quite generally 
in [3]. We recover this result even including the second 
order in the relative temperature difference, which goes 
beyond linear response. For larger deviations from equi- 
librium, however, the efficiency at maximum power in this 
example is smaller than the Curzon-Ahlborn bound. 



^CA = 1 - \/tJ% =%- - 



0{r^) (25) 
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Fig. 3: Efficiency rj* at maximum power as a function of 
(Th — Tc)/Tc for a = 1/2, qca, and 1, corresponding to the 
temperature independent mobility of the example, the Curzon- 
Ahlborn efficiency, and the upper bound of (31), respectively. 
The inset shows the same data in a log-plot. 



General Case. — We now relax the restriction to one 
degree of freedom in a harmonic potential and go back 
to the general case where we also impose two bound- 
ary distributions: Pa(x) at r = (and at t = ti -I- ^3) 
and Pbi'x.) at T ~ ti for the two isothermal transitions. 
As seen from Ec^. (7), minimizing the total mean work 
W (leading to a maximal work output) subject to such 
boundary conditions, which leave the AS term invariant, 
is equivalent to the minimization of the irreversible work 
(3) for both isothermal transitions. Using a normalized 
time f = (t — ti)/{tf — ti) and thus a time-dependent 
distribution 7'(x, f), it follows from (1) that the current 
can be expressed as j = J'(x, f ) /{tf — t./). The irreversible 
work (3) can then be written as 



1 



df 



Minimizing the irreversible work Win- thus is equivalent 
to minimizing the integral in (26) which only depends on 
V{k,t) and ^/(x, f). The time evolution drP = —V • j 
expressed in the time-normalized quantities imposes the 
constraint dfV = —V • J^. Minimizing the integral in (26) 
subject to this constraint leads to the optimal distribution 
7^*(x, f) which does not depend on the total transition 
time. Since this optimization has to be done for each 
of the isothermal transitions separately, we get optimal 
distributions ^'^(x, r/ti) and "P^x, (T-ti)/tz) for the two 
isothermal transitions which both are independent of the 
total transition times tj, j = 1,3. 

The optimization with respect to maximal work output 
is tractable analytically only if we restrict to harmonic 
potentials y(x, r), as shown above in the example. Oth- 
erwise, the resulting Euler-Lagrange equations are non- 
linear partial differential equations and even the numeri- 
cal solution seems to be a difficult task. For a derivation 
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of our central result, however, we do not need an explicit 
solution. It suffices to observe that Eq. (26) implies a re- 
markable universal feature for the optimal driving scheme: 
The (mean) irreversible work in the isothermal transitions 
scales exactly with the inverse transition time 



(27) 



where A^^'^.^ is independent of the transition times tj{j ~ 
1,3) and becomes a functional only of pa{x), P6(x). and 
the mobility matrix fi. As a direct consequence, the work 
output during a complete cycle is given by 



-W^ - (Th - Te)A5 



4(1) 4(3) 
iir irr 



The heat flux from the hotter bath Th is 



(1) 



(1) 



(28) 



(29) 



For negative AS* < and transition times larger than 
the "stall time" t > t^^, where the engine does not perform 
work anymore, the setup performs as a heat engine with 
efficiency 

-W 

(30) 

approaching the Carnot efficiency yyc = 1 ^ Tc/Th for infi- 
nite cycle time t ^ oo. Maximization of P = —W/ {ti +t3) 
with respect to the transition times ii,t3 leads to a uni- 
versal expression for the efficiency at maximum power of 



Vc 



2 - arjc 



where 



7]* (a) 



(31) 
(32) 



(i) 



This expression for the efficiency at maximal power is our 
main result. In the linear response limit, Th — Tc <C Tc, we 
recover rj* '7c /2, consistent with the analysis in [3]. Be- 
yond linear response, one needs the irreversible action A 
which, in principle, can be calculated for any specific sys- 
tem by solving the optimization problem. In general, due 
to the temperature dependence of the mobility fj., the irre- 
versible action A^^l is temperature dependent. An explicit 
statement becomes possible for an isotropic mobility ma- 
trix /J, in Eq. (1). In this case, a depends only on the tem- 
perature dependence of ji via a[^^ /a[^^ = /z(Th)/^(Tc). 
For a temperature independent mobility, i. e. a ~ 1/2, 
we recover the result (23) from the example above which 
thus holds also in higher dimensions. Only for ^ ~ 1/T, 
corresponding to a = acA = 1/(1 + ^/Tc/Th), we recover 
the Curzon-Ahlborn eflfciency rj* = 1 — ^Jt\JT^. These 
efficiencies at maximum power are shown in Fig. 3. From 
a theoretical perspective, it is important to note that for 
1 > a > acA, the result (31) exceeds the Curzon-Ahlborn 



bound. While in the present setup this can be achieved 
only if the mobility [i decreases faster than ^ l/T with 
temperature, a = 1 can be reached in discrete ratchet heat 
engines [25]. 

As outlined in the introduction, the Curzon-Ahlborn 
efficiency is a rather universal bound on the efficiency at 
maximum power for a wide class of model heat engines. 
Why does the efficiency at maximum power derived here 
also under quite general conditions differ significantly from 
the Curzon-Ahlborn result? 

A first hint is given by analyzing our example in anal- 
ogy to an endoreversible Curzon-Ahlborn machine. For 
the latter, a linear Fourier law for the heat flow is as- 
sumed. However, for systems driven strongly out of equi- 
librium, non-linearities can occur. In order to analyze our 
model heat engine in this framework, we first have to as- 
sign an "effective" temperature to the working medium 
(the Brownian particle). Since in this harmonic potential 
the distribution remains equilibrium-like (Gaussian) even 
out of equilibrium, we can assign an effective temperature 



rcff(r) = w(t)\{t) 



(33) 



to the system. If we (instantaneously) coupled a heat bath 
at this temperature Tes to our system and kept the current 
trap strength A constant, we would recover an equilibrium 
situation. We now express the heat exchange of the system 
in terms of this effective temperature 



dQ 



dW 
1 



dE 



1 



-wdX 



1 



diwX) 



1 



-Xdw 



Xwdr = — /iA(T — Xw)dT 
-fiX{T){T - Tcs)dT 



(34) 



and recover a linear Fourier law for the heat flow. 
However, the thermal conductivity /iA(T) becomes time- 
dependent in contrast to the assumptions of the endore- 
versible Curzon-Ahlborn engine. 

Secondly, the derivation of an apparent universal 
Curzon-Ahlborn bound in Refs. [3, 4] by using a cascade 
construction of N machines working in the linear response 
regime comes at the price of an inflnitesimally small power 
output in the (implicitly assumed) limit N —^ oo. For the 
setup in [3], it can easily been shown that for a mechani- 
cal heat engine in the linear response regime, the maximal 
power output scales as P ~ (Ar)^ where AT = Th - Tc. 
The coupling of N such machines with equal temperature 
differences AT = (Th — Tc)/N thus leads to a total maxi- 
mal power output of P* ^ 1/N. 

Role of kinetic energy. — So far, we have neglected 
the kinetic energy of the degrees of freedom. However, 
even in the overdamped limit, coupling to and decoupling 
from the heat baths leads to a heat flux due to the ki- 
netic energy of the particle [13]. This additional heat flux 
Qkin = (Th— Tc)/2 per degree of freedom can easily be con- 
sidered, leading to a general decrease of efficiencies. For 



Tim Schmiedl and Udo Seifert 



large changes of the potential V{x,t) during the cycle, 
however, this contribution is small and can be neglected. 
In general, a given amount of heat leakage Q^'^'^^ per cycle 
(e. g., due to the kinetic energy) leads to an additional 
term in Eq. (29) 



.(1) 



(35) 



For Q''^^'^ > 0, Carnot efficiency cannot be reached any- 
more and for ii oo, oo, we get the quasistatic 
efficiency 



lcak\ 



(36) 



A straightforward calculation shows that the result (31) 
then is generalized to 



ri*{r) 



'7qs 



(37) 



where rjc now is replaced by the quasistatic efficiency rj^g ■ 

Summary. — In summary, we have analyzed the effi- 
ciency at maximal power for a large class of stochastic heat 
engines. We find a universal law (31) differing substan- 
tially from the prevailing Curzon-Ahlborn bound. Even 
though we here derived the result using the concept of 
stochastic thermodynamics, one might expect that under 
appropriate conditions a "thermodynamic limit" exists, 
bridging the gap between the mesoscopic system discussed 
here and macroscopic heat engines. 



[16 

[17: 
[is: 

[19 
[20 

[21 

[22 

[23 

[24 

[25 

[26 

[27 

[28 

[29 



Schmiedl T., Speck T. and Seifert U., J. Stat. Phys., 
128 (2007) 77. 

Feynman R. p., Leighton R. B. and Sands M., The 
Feynman Lectures on Physics (Addison- Wesley, Reading, 
MA) 1966. 

Buttiker M., Z. Phys. B, 68 (1987) 161. 
Landauer R., J. Stat. Phys., 53 (1988) 233. 
Blanter Y. M. and Biittiker M., Phys. Rev. Lett., 81 
(1998) 4040. 

Jarzynski C. and Mazonka O., Phys. Rev. E, 59 (1999) 
6448. 

ASFAW M. and Bekele M., Eur. Phys. J. B, 38 (2004) 
457. 

Gomez-Marin A. and Sancho .1. M., Phys. Rev. E, 74 
(2006) 062102. 

VAN DEN Broeck C, Kawai R. and Meurs p., Phys. 
Rev. Lett, 93 (2004) 090601. 

Velasco S., Roco J. M. M., Medina A. and Hernan- 
dez A. C, J. Phys. D, 34 (2001) 1000. 
Reimann p., Bartussek R., Haussler R. and Hanggi 
P., Phys. Lett. A, 215 (1996) 26. 

SOKOLOV I. M. and Blumen A., J. Phys. A: Math. Gen., 
30 (1997) 3021. 

KULIC I. M., Thaokar R. and Schiessel H., J. Phys.: 
Condens. Matter, 17 (2005) S3965. 

Schmiedl T. and Seifert U., Phys. Rev. Lett, 98 (2007) 
108301. 



REFERENCES 

[1] CuRZON F. L. and Aiilborn B., Am. J. Phys., 43 (1975) 
22. 

[2] Hoffmann K. H., Burzler J., Fischer A., Schaller 
M. and Schubert S., J. Noneq. Thermodyn., 28 (2003) 
233. 

[3] VAN DEN Broeck C, Phys. Rev. Lett, 95 (2005) 190602. 
[4] DE CiSNEROS B. J. and Hernandez A. C, Phys. Rev. 

Lett, 98 (2007) 130602. 
[5] JijLlCHER F., Adjari a. and Frost J., Rev. Mod. Phys., 

69 (1997) 1269. 
[6] ASTUMIAN R. D. and Hanggi P., Physics Today, 55(11) 

(2002) 33. 

[7] Reimann P., Phys. Rep., 361 (2002) 57. 
[8] Sekimoto K., J. Phys. Soc. Jpn., 66 (1997) 1234. 
[9] Sekimoto K., Prog. Theor. Phys. Supp., 130 (1998) 17. 
[10] Sekimoto K., Takagi F. and Hondou T., Phys. Rev. 

E, 62 (2000) 7759. 
[11] Parrondo J. M. R. and Cisneros B. J. D., Applied 

Physics A, 75 (2002) 179. 
[12] Parmeggiani a., JijLlCHER F., AjDARi A. and Frost 

J., Phys. Rev. E, 60 (1999) 2127. 
[13] Hondou T. and Sekimoto K., Phys. Rev. E, 62 (2000) 
6021. 

[14] VAN den Broeck C, Adv. Chem. Phys., 135 (2007) 189. 
[15] Seifert U., Phys. Rev. Lett, 95 (2005) 040602. 



